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Giant magnetoresistance of edge current between fermion and spin topological systems 
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A spin-1 subsystem conjoined along a cut with a subsystem of spinless fermions in the state of 
topological insulator is studied on a honeycomb lattice. The model describes a junction between 
a 2D topological insulator and a 2D spin lattice with direction-dependent exchange interactions in 
topologically trivial and nontrivial phase states. The model Hamiltonian of the complex system 
is solved exactly by reduction to free Majorana fermions in a static Z 2 gauge field. In contrast 
to junctions between topologically trivial phases, the junction is defined by chiral edge states and 
direct interaction between them for topologically nontrivial phases. As a result of the boundary 
interaction between chiral edge modes, the edge junction is defined by the Chern numbers of the 
subsystems: such the gapless edge modes with the same (different) chirality switch on (off) the 
edge current between topological subsystems. The sign of the Chern number of spin subsystem is 
changed in an external magnetic field, thus the electric current strongly depends both on a direction 
and a value of an applied weak magnetic field. We have provided a detailed analysis of the edge 
current and demonstrate how to switch on (off) the electric current in the magnetic field. 


INTRODUCTION 

The Kitaev honeycomb model [l] has given a great 
boost to the investigation of topological states in mag¬ 
netic systems. A gapped phase with a chiral Majorana 
edge state and non-Abelian anyonic excitations in a bulk 
are realized in the model. The ideas stated in ref. |T[ have 
interesting physical applications, ranging from an exact 
method of solution of the models generalized to a vari¬ 
ety of other lattices 0E] to a study of topological chiral 
quantum spin liquids with gapless chiral edges states and 
spin liquids with Fermi surface. In the model of ref. jl] 
the gauge fields are constants of motion, this fact being 
the key to the exact solvability of the model. Physical re¬ 
alizations of the spin- \ Kitaev model have been proposed 
in optical lattice systems [ 4 ]. 

An interest in the topological systems is inspired by 
special properties of their boundaries U®. Excitations 
in the bulk are gapped in topologically trivial and non¬ 
trivial phases (in a weak magnetic field) and at the same 
time boundary excitations are gapped (in topologically 
trivial phase) or chiral gapless (in topologically nontriv¬ 
ial phase). In contrast to the junction between subsys¬ 
tems in topologically trivial phases, the junction between 
topological phases is leaded to a strong interaction be¬ 
tween chiral edge states on the boundary. As a result, 
edge modes are deformed due to the boundary interac¬ 
tion. The unique properties of the surface are responsible 
for their exotic electromagnetic properties, and might be 
used to realize new topological phases hosting Majorana 
modes when brought into contact between different topo¬ 
logical systems. 

The aim of the paper is to demonstrate ‘the work’ of 
the junction between subsystems with different nature 
(spin and fermion subsystems) in topologically nontrivial 
states. We present certain exact results for chiral edges 
modes that define the chiral current along the cut for 
an arbitrary value of a coupling constant. The junction 


considered is unique in the sense that it the first exact 
description of the junction between different subsystems 
in topologically nontrivial states. We propose a princi¬ 
pal new realization of Kitaev’s ideas considering a spin 
sublattice connected by a junction with a topological in¬ 
sulator (TI), in other words, a junction between topolog¬ 
ical spin and fermion insulators. As a model system, the 
topological Kondo lattice thus identifies a new class of 
spin-charge insulator. The experimental situations when 
the spin-charge coherent transport is observable in the 
quantum spin Hall regime of a two-dimensional TI [ 7 ] is 
discussed in ref. [8]. 


MODEL 

We focus on a 2D model of spinless fermions in a topo¬ 
logical state conjoined via a contact interaction with a 
spin-^ subsystem with direction-dependent exchange in¬ 
teractions. The model is defined on a honeycomb lattice 
with a spin- \ located at each site (see fig. Ih). The total 
Hamiltonian Pi = Pif + PL S + Pi cu t describes fermion PLf 
and spin PL S subsystems conjoined via the boundary term 
Picut along the cut. We use the Kitaev model [1] for the 
spin subsystem 

Us=I x X a i a j +I V J2 (T V i <J V j+Iz X a i a h 
x —links y — links z — links 

where i is the nearest neighbor of j connected by x -, y- 
or z-links, crj are the three Pauli operators at site j and 
J 7 is an exchange integral along a link 7 = x, y, z. 

For the sublattice of the spinless fermions we will use 
the following tight-binding model that takes into ac¬ 
count hoppings between the nearest- and next-nearest- 
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FIG. 1. (Color online) Spin and fermion sublattices defined 
on a honeycomb lattice consisting of the spin (blue circles) 
and fermion (red circles) sites, the unit cells highlighted by 
gray box, 60° coordinate system, illustration of the zig-zag 
cut and the interaction term A. 


neighbors with equal hopping and pairing amplitudes 

H ( a \ a m + a l a m ) 

(l,m) 

-it 2 Fi ( a l a m + a l a mJ +h.C., ( 1 ) 

«Z,m» 

where a\ and ai are the spinless creation and annihila¬ 
tion operators defined on the honeycomb fermion sublat¬ 
tice and satisfying the usual anticommutation relations, 
t\ is an overlap integral of the nearest-neighbor hop¬ 
ping between neighbor sites (Z,m) from l of the type A 
to m of the type B (A and B are the atoms of unit cell 
of the honeycomb lattice), £2 is a clockwise (relatively 
honeycomb) the next-nearest-hopping between second- 
neighbors ((Z,m)). 

We have used the notion of Majorana fermions di = 
a t + aj and gi = | (a l — aj ) for the Haldane Hamiltonian 
without a local staggered potential with £2 —^ £2 exp (£</>) 
in the special points (j) = ±7t/2 [10 ] obtaining 

l~Lf — Bf(d) — iti ^ ^ didm ££2 ^ ^ Vimdidmi 

(l,m ) ((Z,ra)) 

where vij = ±1 stands for clockwise (anticlockwise) next- 
near-neighbour hopping inside corresponding plaquette. 
In spite of the availability of the pairing terms in ([!]), the 
Fermi energy of electrons is well defined and the spectrum 
of the Hamiltonian T-L / which is a half of the one-particle 


Haldane Hamiltonian 

Ho = H f (d) + H f ( g ) = -ih ^ d\d m - it 2 ^2 v inA~ a m 

(hi) «bi)) 

which consists of two noninteracting Hf-like parts of d 
and g Majorana fermions, di is a fermion annihilating 
operator. 

The term H cu t is governed by the interaction between 
local moments and itinerant spinless fermions located at 
the cut a (see fig. |T]) 

T~Lcut — A — (^) 

a 

where A is the coupling parameter, a denotes the bound¬ 
ary sites connected by two links (z- and x- links) along 
the cut. We consider only zig-zag cut along the direction 
shown in fig. [l] with black arrow whose naive cartesian 
coordinates are ( 1 / 2 , a/ 3 / 2 ). The Hamiltonian ([2]) de¬ 
fines the junction between subsystems along the zig-zag 
cut with broken ?/-links in sites on the cut. 

The energy of external magnetic field H can be added 
to the total Hamiltonian of the system Hh = — H JT 
where an external magnetic field acts on all spins. The 
system is stable due to a time reversal (TR) symmetry, 
the energy spectrum of spinless fermions and Majorana 
fermions can not open gaps without breaking TR symme¬ 
try. In the case of noninteracting subsystems, an external 
magnetic field in the Kitaev model pQ and hopping be¬ 
tween the next-nearest-neighbor fermions in the Haldane 
model m break TR symmetry of the corresponding sub¬ 
system. The boundary term © breaks TR symmetry 
of the total Hamiltonian F£, but does not stabilize the 
spin subsystem; being the boundary term, it only ‘de- 
forms’ the edge states. Below we consider an effect of 
the boundary interaction on edge states’ structure in the 
charge and spin sectors. 


TOPOLOGICAL EDGE CURRENT 

We follow Kitaev [T using the Majorana fermion repre¬ 
sentation of spin -1 operators within an enlarged Hilbert 
space. The spins are located at the vertices connected by 
three types of links, each link is associated with a type of 
the interaction between nearest-neighbor spins pQ. The 
model is exactly solvable in the absence of an external 
magnetic field, the Hamiltonian can be diagonalized us¬ 
ing a presentation of the Pauli operators crj = ibjcj in 
terms of Majorana fermions b J and Cj satisfying a Clifford 

algebra: {ci,cj} = 2%, {b],^ } = 2< 5 77 «%, {ci,6j} = 0. 

The total Hamiltonian can thus be rewritten in the 
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following form 

A] jCi Cj — it\ did m 

7 =x,y,z < i,j) (l,m) 

-it 2 ^2 v lrn did m + g a d a c a b v a , (3) 

«i,m» « 

where the matrix A consists of A]- = 2I~u]? for the 

O' 

directed links 7 = x,y,z. 

In the absence of an external magnetic field, the op¬ 
erators uJ7 = ib]b^ and wi- = ib^dj commute with the 

LJ L J (JCJ UL J 

Hamiltonian, they are constants of motion with eigen¬ 
values d=l, the model has an exact solution. The state 
with vortex-free field configuration is favored. Numeri¬ 
cal calculations show that the interaction (along the zig¬ 
zag cut) between subsystems does not destroy the uni¬ 
form flux (the ground state) sector. Let us simplify the 
problem and consider the magnetic subsystem with an 
isotropic exchange interaction I x = I y = I z = / = 1. 

Consider the coupled subsystems in the presence of a 
magnetic field counted with perturbation theory accord¬ 
ing to Kitaev’s approach, an effective magnetic field h 
is defined as H x H y H z /1 2 - H 3 /1 2 jj. We solve eigen¬ 
value problem numerically and calculate energies of ex¬ 
citations in the sectors with uniform flux. The fermion 
subsystem is TI for arbitrary values of £2 mm, whereas 
the ground-state phase diagram of the spin system is de¬ 
fined by the values of exchange integrals. In the case 
of an isotropic exchange interaction the phase state of 
the Kitaev model [L is topologically nontrivial one with 
Chern number 1. Spin and fermion excitations with chi¬ 
ral edge modes for the zig-zag cut are shown in fig. [2] for 
different topological states of spin subsystem. The sign 
of the Chern number of TI varies with the sign of £2 m, 
whereas the sign of the Chern number of spin subsytem 
is defined by the sign of an external magnetic field pQ. 

The structure of edge modes depends on the Chern 
numbers of noninteracting subsystems: due to the 
boundary interaction the chiral gapless edge modes are 
‘deformed’ into gaped modes in the case of different 
Chern numbers (see fig. §b)) and into gapless modes 
for the same Chern numbers (see fig. [2|d)). In contrast 
to ref. [7], the chirality of the edge modes is conserved. 
The chiral current along the cut is equal to zero (cur¬ 
rentless state) in the first case, a dissipativeness current 
along the cut (current state) is realized in the second 
case. The Chern number of the spin subsystem depends 
on the direction of a magnetic field, therefore the current 
state of the junction is defined by an applied magnetic 
field. In order to change the value of the current along 
the cut between subsystems, we can change the direction 
of the magnetic field, in other words, an electric current 
along the cut depends anomalously on an applied mag¬ 
netic field. 




FIG. 2. (Color online) Dispersion s(k) of the system cal¬ 
culated on two cylinders with open boundary conditions, a 
junction for zig-zag boundaries with t\ = 5, £2 = 1, / = 1, 
A = 0.5 a) e{k) of bulk and edge modes along the zig-zag 
direction for h = 0.2 b) the edge modes for h = 0.2, c) e(k) 
for h — —0.2, d) the low energy energy levels h — —0.2. Pure 
fermionic states are marked with red color, pure spin states 
are marked with black, mixed states are shown by gradient 
through orange (fermion part dominating) to blue (spin part 
dominating). 

Due to the hybridization of the edge modes, pure 
charge and spin chiral currents are not realized on the 
boundary, spin and charge currents are strongly corre¬ 
lated. We have calculated partial (spin and charge) cur¬ 
rents as functions of the parameter of the interaction and 
the value of magnetic field, and displayed them in fig. [3] 
The charge edge mode is high-energetic one, therefore the 
spin current dominates in low energy region (see fig. [ 3 ]) . A 
main result is presented in fig.^a) and fig.^b) where the 
charge current jp and magnetoresistance R = (^^) _1 
calculated as a function of magnetic field, fig. [3ji and 
fig. J^b) illustrate the results obtained before in fig. [ 2 ] 
and show an anomalous dependence of chiral edge charge 
current on value and a sing of magnetic field. Low- 
temperature magnetoresistance of one-dimensional sur¬ 
face transport has been observed in topological Kondo 
insulator SmB 6 [9j. 

CONCLUSIONS 

We have considered a 2D complex topological sys¬ 
tem — spin subsystem (in the framework of the Kitaev 
model) and TI (in the framework of the Haldane model 
in special points without a staggered potential) which 
are in contact along the boundary and calculated the 
chiral current along the boundary. We emphasize the 
fact that results obtained are exact, the models proposed 
have exact solution, therefore all the kinetic processes 
are dissipativeness. We have shown that the interaction 
on the cut leads to a strong hybridization of spin and 
charge edge modes and changes cardinally the transport 
of spin and charge in gapless edge states. We have cal- 
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FIG. 3. (Color online) Fermionic (red curve) and spin (black 
one) currents along the cut and their ratio depending on 
model parameters A and H for ti = 5, 1=1, £2 = 1 : a) 
currents for A = 0.2 with H varying, b) derivative R~ x = 
depending on H, c) currents ratio depending on H for A = 0.2, 
d) currents ratio depending on A for h = 0.2. 


ciliated the structure of chiral edge modes in topolog¬ 
ically nontrivial phases of subsystems and partial edge 
currents. Due to the hybridization on the boundary, the 
gapless edge modes with different chirality are attracted 
and connected to one gapped edge mode. At the same 
time the gapless edge modes with the same chirality are 
repulsed conserving gapless states. The state along the 
boundary is currentless in the first case and there is a 
current in the second case, the transition from one to an¬ 
other state is realized in weak magnetic field by changing 
a direction of an applied magnetic field. We have shown 
that topological properties of interacting systems with 
different nature (charge and spin) can define anomalous 
magnetic resistivity chiral edge current and lies to giant 
magnetoresistance effect. It would be of interest to find 
materials which exhibit this effect, as well as possible 3D 
generalizations. 
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